Anomalous Conductance Suppression in Magnetic Clusters 
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Recent scanning tunneling microscope measurements on magnetic clusters reveal an anomalous 
conductance dip at zero bias voltage, whenever the cluster consists of an odd number of magnetic 
atoms [Science, 312, 1021 (2006)]. We address the conductance anomaly within a model, and 
reproduce the conductance dip with the only assumption that the magnetic cluster has a conducting 
level near the Fermi level of the system. We show that the width of the conductance dip scales 
quadratically with external magnetic fields, in excellent agreement with experiments. In addition, 
we predict that the presence of the impurity will be measurable as inelastic Friedel oscillation in the 
substrate density of electron states. 

PACS numbers: 72.25.-b, 73.23.-b, 73.63.Kv 



The study of spin-dependent transport phenomena in 
mesoscopic systems is becoming a well established re- 
search field [l| which has gained lots of attention recently, 
both experimentally Q and theoretically While one 
major reason is the prospect of using such systems in 
spin-qubit readout , fundamental questions concerning 
how spin coherence and interactions influence the trans- 
port are of main interest. Measurements of the effects can 
be performed by, for example, electron spin resonance 
(ESR) [6[, and more recently, ESR-scanning tunneling 
microscopy (ESR-STM) techniques 0,0. In this context, 
the ESR-STM is especially interesting because of the pos- 
sibility to manipulate single spins [3 HI U | , something 
which is crucial in spintronics and quantum information. 
Experimentally, modulation in the tunneling current has 
been observed by STM using spin-unpolarized electron 
beam 0, 0] ■ Lately, there has also been a growing inter- 
est in using spin-polarized electron beam for direct de- 
tection of spin structures [l2| , as well as utilizing the in- 
elastic electron scanning tunneling spectroscopy (IETS) 
for detection of local spatial variations in electron-boson 
coupling in molecular systems 0, [wj- IETS was re- 
cently suggested to be used for surface imaging of in- 
elastic Friedel oscillations caused by inelastic scattering 
effects on a local impurity [l5j |. 



Recent STM measurements of magnetic atoms, e.g. 
Mn, located on a metallic surface revealed conductance 
dips at vanishing bias voltage whenever the mag netic 
cluster consisted of an odd number of atoms [lg, Il7 |. 
The conductance dips were assigned to inelastic scatter- 
ing, e.g. spin-flip, events between the spin channels in the 
system. In case of an even number of magnetic atoms in 
the cluster the zero bias dip were absent, however, the 
conductance were suppressed for a larger range of bi- 
ases between. The limits of this range were determined 
by the energy barrier needed for spin-flip transitions be- 
tween the singlet (S — 0, m z — 0) ground state and the 
triplet states (S = 1, m z — 0, ±1). 



The central question we address in this paper is the 
anomalous feature at zero bias voltage seen in the STM 
conductance for an odd number of magnetic atoms in 
the cluster. In this case the ground state of the mag netic 
cluster, which is anti-ferromagnetically ordered |l8|, is a 
doubly degenerate spin state. Due to (inelastic) spin-flip 
transitions in the tunneling between the magnetic cluster 
and the STM tip and the substrate, further channels open 
for conduction in addition to the existing spin-preserving 
channels. The result is that the tunneling electrons be- 
come shared between different pathways in the system, 
see Fig. Q]a), which gives rise to Fano-like interference 
effects between tunneling electron wave functions. The 
interference tends to suppress the conductance through 
one of the channels whereas the conductance in the other 
is enhanced. The destructive interference in the for- 
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FIG. 1: (Color online) a) Sketch of the branching of the tun- 
neling electron wave between the two spin-channels due to the 
strong spin-flip tunneling rate, b) Quadratic fit of the FWHM 
of the experimental conductance dips, extracted from Ref. 
[l7j . as function of the external magnetic field, c) Calculated 
differential conductance for different external magnetic fields 
within the present theory, e.g. using Eq. (0 with equal spin- 
preserving (Fq) and spin-flip (r s ) tunneling rates (r s = To). 
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mer case, generates a markedly suppressed conductance 
around zero biases whenever the spin-flip tunneling rate 
is of the order of the spin-preserving tunneling rate. 

The result of this Letter is that we reproduce the con- 
ductance dip within a model with the only assumption 
that the magnetic cluster has a conducting level at the 
Fermi level of the system. Through the model, we find 
that the width of the conductance dip scales quadrati- 
cally with external magnetic field B, in excellent agree- 
ment with experiments, c.f. Fig. [T]b). 

We consider a model where electrons tunnel between 
the tip and the substrate via a molecular level e a , as- 
sumed to be close to the common equilibrium chemical 
potential Ef of the tip — atomic cluster — substrate- 
system. We argue that we can easily reproduce the dip 
in the tunneling conductance for odd spin state within 
this model, with the only assumption that e a < Ef- For 
simplicity, we assume that there are no electron-electron 
or spin-flip interactions in the metallic leads. The spin- 
flip transitions are confined in both barriers between the 
local spin and the leads. The system is modeled by the 
following Hamiltonian: 



U = ^2 e ka cl a c ka + ^ Zod\d G 

kcr <y 

+ ^(<w4A' + #-c-), (i) 

where ct (cka ) creates (annihilates) an electron with mo- 
mentum k, spin a =f , |, and energy £ k(T in the left (L) 
or right (R) lead, whereas d\ (d a ) creates (annihilates) 
an electron with spin a at the energy e a = e + crA/2 
on the local spin site, where A is the spin split on the 
impurity spin, which may be caused by external mag- 
netic field B and/or the internal structure of the mag- 
netic cluster. Electrons on the spin site hybridizes with 
those in the leads with the rate v kaa i ■ Consider the 
case of no spin-orbit interaction but purely exchange 
interacting atomic orbitals with single electronic spin, 
e.g. JS a t om ic • Scicctronic- The Fermi golden rule already 
would indicate that the lifetime broadening r CT of such 
electronics states as well as the spin-flip scattering rate 
T s all are proportional to J 2 [Nf(ef) + Nr(ef)], where 
Nli r {ef) is the Fermi level density of states (DOS) in 
the tip/substrate. 

The current through the system is derived by stan- 
dard methods, considering the temporal derivative of the 
occupation Nl = ^2 paeL cp a c pcr in the tip (left lead), 
i.e. I = ~ed/dt(NL}- The differential conductance 
G{V) = dl/dV can then be written as [3] 



em = e -£ «r 



T(u) cosh" 2 0(cj - eV)duj, (2) 



where T{lu) = T L G r (uj)T R G a (uj) is the transmission 
coefficient. Here, the coupling between the local spin 
and the tip/substrate is given by T L I R , with T^J 1 = 



27r J2ka"eL/R v ka"a v kcr"a'S{uJ-Ekcr") 1 whereas G r/a {uj) is 
the Fourier transformed retarded/advanced Green func- 
tion (GF) for electrons on the local spin site: G r J^,(t) = 
Ti0{±t){{d a {t),dl,(0)}). Note that both T L ' R and G r ' a 
are matrices in the spin space of the impurity. 

The impurity site GF is solved by equation of motion, 
which yields (defining A a = A aa for shorter notation) 



GUu>) 



(cu — uj+)(lo — UJ T _) ' 

V'' 

(lo — o;!^)(cij — uj r _ ) ' 



(3a) 
(3b) 



where the self-energy matrix due to the couplings to the 
leads is 



yr _ 



\p V ka n a Vko"a' 

^ W - Eka" + M + ' 
ka"GLL)R 



(4) 



The GFs and G r ^ are obtained by exchanging 
in the above equations. For simplicity we assume large 
band-witdths in the leads, hence, we can approximate 
the self-energies by ££ CT , = -in J2ka" v ka" <y v k°" - 
£ka") = (^"i/^)(^aa' + ^aa 1 )- ^ n tne present case, where 
the magnetic atom is located on a substrate, the cou- 
pling between the impurity site and the tip (left lead) 
is much smaller than the coupling between the substrate 
and the impurity site, i.e. T^ a , /T R a , =A<Cl. Hence, 
the self-energies for the impurity site electrons are then 
dominated by the coupling to the substrate (right lead), 
hence, S^, « — iT R a ,/2. Hermiticity of the couplings 

also require that T+( R — — T^/ R . In the experi- 

mental setup, the tunneling current is spin-independent 

[ana. 

Hence, it is reasonable to let T L = Tq and 
T R = T R , where = XT and T R = T . In this case, 
the poles uj± of the GFs in Eq. ([3]) are given by 



»(r t VtI^a^/2. 



(5) 



Suppose that the spin-split A = g^isB, where g and 
fj,B is the gyromagnetic ratio and Bohr magneton, respec- 
tively. We assume that A <C ro,r,,, which is reasonable 
for tunneling rates of the order 10 meV. In this case, the 
widths of the poles scale quadratically with the magnetic 
field, since the width can be written 



r T-r s 



r T-r s 



(G) 



As seen, in the absence of magnetic field B, the pole oj + 
acquires a vanishing width as T s — > To, while the width 
of the pole w_ appoaches To, which is also clear from Eq. 
([5]) . Physically this means that an electron in the level uj+ 
only weakly interacts with the conduction electrons in the 
leads. Therefore, the conduction through the magnetic 
cluster tends to decrease in the limit r, — *■ T,, since then 
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the electron occupation of the level u>+ approaches unity 
whenever the localized level £o < £f-- 

Hence, the conditions that yield the vanishingly small 
width of the pole uj + give rise to an anti-resonance in 
the transmission T{uS) at to = sq. This is easily seen by 
substituting the expression for the GF G r / a {uj) on the 
local spin site into the expression for the transmission, 
which gives 



A 

4(w - 



r 2 
1 



r 2 Q 



•r 



■(A 2 



r 2 

1 



(7) 



When the spin-flip tunneling rate is comparable to the 
spin- preserving tunneling rate, e.g. T s « To, the last 
term in this expression vanishes, leading to the transmis- 
sion 



3» 



(c-e ) 4 + (c- £ 2 )(r 2 -A 2 /2) + (A/2r 



(8) 

which exactly vanishes for uj = Eq and all finite magnetic 
fields. Thus, assuming that the local level £o — 0, the 
differential conductance of the systems will display a dip 
for vanishing bias voltages, see Fig. Q]c). Absence of 
magnetic fields leads to A = 0, which reduces Eq. ([5]) to 



1 



(^-e ) 2 + r 2 ' 



that is, a single peak around Sq of width Tq. The width 
of the anti-resonance in the transmission is related to 
the width of the poles lo+. Thus, the width of the anti- 
resonance scales quadratically with the magnetic field ac- 
cording to Eq. ([6]), something which is also seen in Fig. 
[T]c). This feature is in excellent agreement with the ex- 
perimental findings in Ref. [ItJ ■ 

The dip in the transmission can be understood as an ef- 
fect of destructive interference of the branched tunneling 
electron wave functions between the two spin channels in 
the system. This interference effect also becomes visible 
in the local DOS in the substrate caused by the presence 
of the localized magnetic atom. In the absence of the 
magnetic impurity the local DOS of the substrate is fea- 
tureless. In terms of GFs, the unperturbed DOS is given 
byN R = -(l/vr)Im£ ?CT 3 g r CT M = J2 qa 5(u-e q<r ), where 
g r qa {u>) is the Fourier transform of the retarded substrate 
GF g qa {t) = (-^(Tcg^c^O)} is the retarded unper- 
turbed GF of the substrate. The presence of an impurity 
on the substrate, will then manifest itself in the local 
substrate DOS as a mirror image of the local DOS of 
the impurity, since the substrate GF, now redefined as 
9qaq'a'{t) = (-£) (Tc g(T (t)c\, „, (0)) , becomes dressed by 
the presence of the impurity. Accordingly, the Fourier 
transformed retarded substrate GF can be factored as 



dqaq'a' ~ &qq'&oo' 9qcr + 9 qa v qaa\ ^ aia 2 V q' a' a 2 9q 



(9) 



where we assume summation over repeated indices. The 
total DOS is then given by N r (oj) = N R + SN r {lu), 
where 



8N R (w) w -tr lmT R G r (uj) 



(10) 



= (r 2 + r 2 s )(u - £o ) 2 + |(rj - r 2 )(r 2 - r 2 + a 2 ) _ 

2\(tU — LU + )(u> — LU-)\ 2 

Furthermore, in the limit T s = Tq we find that the fre 
quency derivative d u 8N (u>) is given by 

d 



du> 



5N R (cu) = -2r 2 



(11) 



(c- £o )[(c- £o ) 4 + (A/2) 4 



[(u - e ) 4 + (w - e ) 2 (r 2 - A 2 /2) + (A/2) 4 ] 2 ' 

Clearly, the presence of the impurity provides a finger- 
print in the local substrate DOS, such that the density 
of the impurity manifests itself as a dip around lu = eo in 
the substrate DOS as shown in Fig. [2] a) and b). Addi- 
tionally, the anomalous conductance dip will be identified 
as a peak in the substrate DOS exactly at u> — eo, Fig. [5] 
a). This is understood since an electron in the level u>+ 
interacts only weakly with the electrons in the substrate 
and will, hence, have a negligible influence on the sub- 
strate DOS. Therefore, the substrate DOS returns to its 
undisturbed value at u> = Sq. 

The result for the substrate DOS indicates measurabil- 
ity of the conductance dip as inelastic Friedel oscillations 
in real space on the surface. In order to see this, consider 
the real space GF of the substrate using Eq. ((9| 

g(r, r'; w) - ,g (r, r'; ^)+ ff0 (r, 0; W )S( W ) ffo (0, r'; w). (12) 

In two dimensions, the bare GF becomes 

5 o(r,r';c) = ^ ffo (q,c i ;) e ^( r - r ') 
q 

uj + D 



2TT — J Q {k F \Y-Y'\)[\0g 



uj-D 



iir , (13) 



where k F is the Fermi wave vector and Jo (a;) is the zero 
Bessel function, whereas the self-energy is identified by 
E(w) = (l/2)tr T r G t (uj). The correction to the local 
substrate DOS is given by 

SN(t,u) = --Im{ 5o (r,0; W )E( W )g (0,r; U )}. (14) 

7T 

The effects from the inelastic scattering responsible for 
the conductance dip are included in the self-energy, thus 
we evaluate d8N(r,uj)/du> which is directly proportional 
to cPl/dV 2 in the transport measurements. The ex- 
pected sharp features at e is connected to the spin-flip 
tunneling. Noting that Re{go} ~ 0, our calculations sim- 
plify to SN(r,u) « -(l/7r) 5o (r,0;a;)[ImS(w)] ff o(0,r;w), 
which gives, c.f. Eq. (|10p . 



SN{r,u) 



7T 



2m 



J^{k F r)trImr R G r (uj), (15) 
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FIG. 2: (Color online) a) Correction 8N(ro,w) to the local 
substrate DOS, and b) frequency derivative d u 8N(ro, ui) at 
the impurity. Spatial and frequency dependence of c) the 
correction 5N{r,u), and the derivative d u 6N(r,u) at B — 
4 T. The bold lines display the spatial Friedel oscillations 
that would be measured in the substrate DOS, whereas the 
thin lines show duj8N(ri,u}) and du,SN(ri,bj), that would be 
measured at position n. 



where r — |r|, and a corresponding expression for the 
frequency derivative d^SNfc, u>). This shows that the in- 
elastic (spin-flip) tunneling between the substrate and 
the magnetic impurity will be manifested as inelastic 
Friedel oscillations in the spatially resolved substrate 
DOS at the energy w = £q, see Fig. [2]c) and d). 

In summary, within a model we reproduce the anoma- 
lous conductance dips in the STM measurements of mag- 
netic clusters with an odd number of magnetic atoms. 
We show that the conductance dip is caused by a lo- 
cal level in the magnetic cluster which has a unit oc- 
cupation probability and which hybridizes only weakly 
with the conduction bands in the tip and substrate. The 
weak hybridization between the level and the conduction 
bands is generated by strong spin-flip tunneling rate. We 
demonstrate that the isolated local level will manifest its 
presence through inelastic Friedel oscillations in spatially 
resolved STM measurements of the substrate DOS. 
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